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Abstract 

We describe several studies to measure the charged track reconstruction efficiency and asymmetry of the BaBar de- 
tector. The first two studies measure the tracking efficiency of a charged particle using r and initial state radiation 
decays. The third uses the r decays to study the asymmetry in tracking, the fourth measures the tracking efficiency 
for low momentum tracks, and the last measures the reconstruction efficiency of particles. The first section also 
examines ffie stability of the measurements vs BaR4R running periods. 
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1. Introduction 

The BaBar experiment operated from 1999 to 2008 
at the PEP-II asymmetric e^e" collider at the SLAC 
National Accelerator Laboratory. _BaR4R was designed 
to study CP violation and other rare decays in flavor 
physics from events produced at or near the T reso- 
nances, from 9.46 GeV to over 1 1 GeV. A critical re- 
quirement for meeting BaR4R's science goals was the 
ability to efficiently and accurately detect stable charged 
particles, or tracks, produced in e^e^ collisions. Many 
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analyses performed at BaBar require a precise estimate 
of the track finding efficiency, as input for measuring 
the absolute or relative rate of the physics process being 
studied. 

In this paper, we present the algorithms and methods 
used in BaB/^ to estimate the track finding efficiency. 
To cover the range of particle momenta and produc- 
tion environments relevant to most BaHm. analyses, a 
number of methods are used. To compute the track- 
ing efficiency from data alone, these methods rely on 
special data samples, where additional constraints can 
be applied. The primary efficiency result is computed 
using e^e^ — » r^r" events, which can be cleanly iso- 
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24 lated in the BaBar data sample, and which have a sim- 

25 pie topology. To cross-check this result, we indepen- 

26 dently measure the tracking efficiency using radiative 

27 e^e^ n^n'n^ji'jisR events, where jisr is an initial 

28 state radiation (ISR) photon, which can be constrained 

29 kinematically. To study the reconstruction efficiency of 

30 low momentum tracks, we use D*^ — > D^tt* decays. We 

31 also present a dedicated study of the efficiency to recon- 

32 struct K^^ — > n^n' , whose daughter tracks can have a 

33 different efficiency due to their displacement from the 

34 primary event origin. 

35 The strategy for the r-based and e+e" — > 

36 TT^TT^n^n^yisR track reconstruction efficiency measure- 

37 ments is to use charge conservation and kinematics to 

38 deduce the existence of a track, given a subset of de- 

39 tected tracks in well-defined events. The efficiency 

40 analyses based on D'^ decays and for the efficiency 

41 study use a statistical approach, using properties of 

42 momentum distributions which will be described be- 

43 low. Systematic errors are estimated using internal self- 

44 consistency measures and by comparing different effi- 

45 ciency analysis techniques. 

46 The BaBar detector geometry, material, and sensor 

47 response functions have been accurately modeled in a 

48 detailed simulation based on the Geant4 framework. 

49 The output of the BaBar simulation is processed using 

50 the same reconstruction algorithms as applied to data, 

51 and the results have been found to be very similar to 

52 what we see in data. By using accurate computer mod- 

53 els of the physics processes relevant at BaBar energies 

54 10, [stl; we are able to generate equivalent samples of 

55 simulated data as used in nearly all BaBar analyses, 

56 including the tracking efficiency analysis. BaBar has 

57 therefore adopted the strategy of estimating the track- 

58 ing efficiency relative to that observed in the simula- 

59 tion, which simplifies the application of the tracking 

60 efficiency results in analysis. As will be shown in the 

61 following sections, for most of the studies, the tracking 

62 efficiency found in data agrees within errors with the 

63 efficiency found in simulated data. This allows the re- 

64 suit of the tracking efficiency measurement to be used 

65 in analysis simply by propagating the appropriate sys- 

66 tematic errors on the tracks involved to the simulation 

67 estimate of the analysis signal reconstruction efficiency. 

68 This strategy has been used in many scientific BaBar 

69 publications. However, for analyses involving a K^, a 

70 correction is required in the MC for its daughter recon- 

71 struction efficiency which depends on the kinematics of 

72 the decay of interest. 



73 2. BaBar Detector and Data Sample 

74 The BaBar detector is a multi-purpose device de- 

75 signed to simultaneously measure many properties of 

76 the multiple particles produced in e^e^ collisions near 

77 the T resonances, as described in detail in |4l]. Charged 

78 particles are identified in a Silicon Vertex Tracker 

79 (SVT), and a Drift Chamber (DCH), which are sur- 

80 rounded by a superconducting solenoid that generates 

81 an approximately uniform 1.5 Tesla magnetic field in- 

82 side the sensitive volumes of these detectors. 

83 The SVT consists of five layers of double-sided sili- 

84 con strip detectors covering the full azimuthal range and 

85 the laboratory frame polar angle (diai,) range 20.1° < 

86 < 150.2° |i4|] . The intrinsic resolution of individual 

87 SVT position measurements of particles which traverse 

88 it varies between lOyum and SOyum, depending on the 

89 incident particle direction and the readout view. The 

90 DCH consists of 7,104 hexagonal drift cells, which are 

91 approximately 1.9 cm- wide and 1.2 cm-high, made up 

92 of one sense wire surrounded by six field wires. The 

93 sense wires are lOfim gold-plated tungsten-rhenium, 

94 the field wires are 120/im and 80 /im gold-plated alu- 

95 minium. The cells are arranged in 40 cylindrical lay- 

96 ers and the layers are grouped by four into ten super- 

97 layers extending from roughly 25 cm to 80 cm in the 

98 transverse direction, with full coverage over the range 

99 24.8° < 9iah < 141.4°, and partial coverage over the 

100 range 17.2° < Qiab < 152.6°. The intrinsic resolution of 

101 individual measurements of track position in the DCH 

102 varies between 100 /im and 200 //m, depending on the 

103 track position and angle relative to the wire, with an av- 

104 erage resolution of 150 /im. 

105 The BaBar detector includes a dedicated charged par- 

106 tide identification (PID) device based on detection of 

107 internally reflected Cherenkov radiation (DIRC), and 

108 a Cesium-iodide crystal electromagnetic calorimeter 

109 (EMC) for identifying electrons and photons. The steel 

110 for the solenoid magnet flux return is instrumented with 

111 position-sensitive chambers, which produce distinctive 

112 signatures from passing muons and pions. _BaR4R es- 

113 timates the species of charged particles using a com- 

114 bination of information from these devices, plus the 

115 specific ionization (dE/djc) measured in both the SVT 

116 and DCH. By studying the response of these systems 

117 to high-purity control samples, likelihood functions de- 

118 scribing a track's consistency with each of the 5 charged 

119 particle species (e=^,ju'', n"^, K^, and p"^) directly observ- 

120 able in the BaBar tracking system are defined. Samples 

121 of specific particle species of varying efficiency and pu- 

122 rity are selected by cutting on appropriate likelihood ra- 

123 tios. 



124 The results presented in this paper are based on the 

125 full BaBar data sample, collected in seven distinct pe- 

126 riods. Runs 1-7. Runs 1-6 correspond to data collected 

127 with a center-of-mass (CM) collision energy near or at 

128 the T (4S) resonance and Run 7 corresponds to the data 

129 collected with a CM collision energy at the T(35) and 

130 T (25 ) resonances. 

131 3. BaR4R Track Reconstruction Algorithms 

132 Tracks are reconstructed in BaBar using a combi- 

133 nation of several algorithms. Tracks with transverse 

134 momentum above roughly 150 MeV/c are principally 

135 found in the DCH. Track segments are identified as con- 

136 tiguous sets of hits in a super-layer having a pattern con- 

137 sistent with coming from a roughly radial track. Seg- 

138 ments are linked using their position and angle to form 

139 a track candidate. Track candidates are fit to a helix, 

140 which is used to resolve the left-right ambiguity, and to 

141 remove outlier hits. The candidate is kept if at least 20 

142 DCH hits remain. Tracks with large impact parameter 

143 are found in the DCH using a less restrictive algorithm. 

144 Tracks found in the DCH are fit using a Kalman filter 

145 istl fit, which accounts for material effects and corrects 

146 for magnetic field inhomogeneities. The Kalman filter 

147 track fit is extrapolated inwards, and SVT hits consis- 

148 tent with the extrapolated track position and covariance 

149 are added. 

150 Tracks with low transverse momentum are found 

151 principally in the SVT using hits not already associated 

152 with tracks found in the DCH. Sets of four or more (p 

153 hits (which measure the position in the plane transverse 

154 to the beam direction), in different layers of the SVT 

155 and consistent with lying on a circle, are selected. Hits 

156 in the orthogonal (z) view of the same wafers as the (p 

157 hits are then added to form three-dimensional track can- 

158 didates. Candidates with at least 8 hits are selected, and 

159 fit using the Kalman filter. Additional tracks are found 

160 in the SVT using space points constructed from pairs of 

161 <p and z hits not already used in other tracks. Sets of 

162 at least 4 space points consistent with a helix fit are se- 

163 lected as tracks. DCH hits are added to tracks found in 

164 the SVT in a procedure analogous to how SVT hits are 

165 added to tracks found in the DCH. 

166 After all the tracks in an event are found, they are fil- 

167 tered to remove duplicate tracks due to hard scattering 

168 in the material separating the SVT and the DCH, de- 

169 cays in flight, or pattern recognition errors in the DCH, 

170 where stereo and axial hits generated by a single particle 

171 are sometimes reconstructed as separate tracks. A final 

172 pass to remove inconsistent hits and to add individual 



173 hits missed in the pattern recognition is then performed 

174 using the Kalman filter fit. 

175 The resultant set of tracks is referred to as Charged 

176 Tracks (CT). A Good Tracks (GT) subset of tracks, with 

177 a higher probability of originating from the primary 

178 e^e" interaction, is selected from these. The GT se- 

179 lection requires the impact parameter with respect to 

180 the average interaction point be less than 1.5 cm in the 

181 transverse direction, and less than 2.5 cm along the ma- 

182 gentic field (z) direction. Analyses at J3aR4R generally 

183 use either the CT or the GT track selection, and the 

184 tracking efficiency studies described in this note are per- 

185 formed independently for both. 

186 4. Tau31 Tracking Efficiency Study 

187 The efficiency of charged track reconstruction at 

188 BaBar is determined using e^e^ — > t^t" events. With 

189 over 430 million t pair events collected at BaBar, t de- 

190 cays provide an opportunity to make a precision mea- 

191 surement of the tracking efficiency. At the CM energies 

192 produced at BaBar, t decays are an ideal candidate for 

193 measuring the tracking efficiency because they have a 

194 momentum and angular distributions of tracks that are 

195 similar to those from decays of D and B mesons. De- 

196 cays of T leptons have a high track density due to the 

197 initial boost, fi ~ 0.94c, of the r leptons, while the total 

198 track multiplicity is low. The r lepton has a life-time 

199 of (290.6 ± 1.0) X 10^'^i, which results in a transverse 

200 flight length of 200 jum at the BaBar CM energies, a 

201 value that is slightly larger than the beam spot size but 

202 small enough not to impact the tracking efficiency. 

203 The tracking efficiency is measured using e^e^ — » 

204 T^T" events in which one r lepton decays leptoni- 

205 cally via — > yu^v^v^, and the other t lepton de- 

206 cays semi-leptonically to 3 charged hadrons via r* — > 

207 h^li^h^Vj + > neutrals (excluding K'^), referred to 

208 as Tau31 events. The tracking efficiency is measured 

209 using the 3-prong t decays. The branching ratio of 

210 — > is^VjjVt and 3-prong t decays are (17.36 + 0.05)% 

211 and (14.56 ± 0.08)% M respectively, so that Tau31 

212 events constitute over 5% of the total. The t pair candi- 

213 dates are selected by requiring an isolated muon track, 

214 plus at least two other reconstructed tracks consistent 

215 with being hadrons. Events are selected in two over- 

216 lapping channels; those where two of the hadrons have 

217 the same charge ("same-sign"), and those where two of 

218 the hadrons have opposite charge ("opposite-sign"). Re- 

219 quiring a muon track is an essential part of suppressing 

220 non-T backgrounds: radiative Bhabha events where the 

221 photon interacts with the detector material producing 

222 an e^e^ pair (conversion), y-y events, and qq events. 
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Figure 1 : The trae cos(6) in the laboratory frame for the fourth track 
for the selected opposite-sign and same-sign MC events. N = N3 + 
N4 + N5 is the number of selected same-sign and opposite-sign events. 
A'4 is the number of events where the fourth track is found for the 
CT definition, N3 is the number of events where the fourth track is 
not found for the CT definition and is the number of events where 
two CT candidates are found for the fourth track. The dotted fines 
indicate the outer edge of the tracking detectors, while the dashed 
lines indicate the edge where the full detector coverage begins. The 
region in between where there is partial coverage is indicated by the 
shading. 



223 Charge conservation infers the existence of the fourth 

224 track. 

225 The tracking efficiency e is defined by 
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226 where A is the geometric acceptance of the fourth track 

227 constrained by the t pair kinematics and the selection 2=3 

228 criteria of the Tau3 1 sample, A^^ is the number of events 254 

229 where the fourth track is found, and A^3 is the number 255 

230 of events where the fourth track is not found. The ge- ^'^ 

231 ometric acceptance of the BaBar detector for a uniform 

232 cos{9) distribution is ~ 83%. In figures [T] and |2] the 2=" 

233 geometric acceptance of the detector is plotted for sim- 259 

234 ulated events as a function of the polar angle (6) and the 260 

235 transverse momentum (p,) of the fourth track, respec- 261 

236 tively. These figures demonstrate the limited angular 262 

237 acceptance of the detector, and the poor acceptance for 253 

238 low momentum tracks. 254 

265 

239 4.1. Monte Carlo Samples 

267 

240 r^r" pair events are simulated with higher-order ra- 268 

241 diative corrections using the KK2f Monte Carlo (MC) 259 

242 generator |01 with r decays simulated with Tauola fsl, 270 

243 |9[]. The simulated Standard Model backgrounds in- 271 

244 elude: bb; cc; ss; uu\ and /i^ju" events 272 



Figure 2: The true p, for the fourth track for the selected opposite- 
sign and same- sign MC events. N = N3 + N4 + is the number 
of selected same-sign and opposite-sign events, ^4 is the number of 
events where the fourth track is found for the CT definition, ^3 is 
the number of events where the fourth track is not found for the CT 
definition and Nf is the number of events where two CT candidates 
are found for the fourth track. The tracks in the shaded region do not 
reach the outer edge of the DCH. 



The number of simulated background events is compa- 
rable to the number expected in the data, with the ex- 
ception of Bhabha and two-photon events, which are not 
simulated. Bhabha and two-photon events backgrounds 
are studied with control samples. The detector simula- 
tion and reconstruction of the MC events is described in 
Section|2] 

4.2. Event Selection 

We require the events to have a minimum of three GT 
and a maximum of five CT tracks. Events with K^^ are 
removed, where the candidate is defined as having 
two oppositely charged tracks with an invariant mass 
within 10 MeV of the mass a vertex displaced 
more than 2 mm from the beam-spot and a vertex fit x' 
probability of more than 1%. The three GT tracks are 
required to have p, > 100 MeV. To remove any remain- 
ing duplicate tracks, the three GT tracks are required to 
satisfy an isolation cut in 9, (p and momentum by 0. 1 
rad, 0. 1 rad and 0.4 GeV, respectively. One of the three 
GT tracks must be more than 120 degrees from the other 
track. This isolated track must satisfy a tight muon PID 
selection. At least two of the other tracks are required 
to be identified as pions, by being inconsistent with a 
loose electron PID selection. 

For the "same-sign" channel (t'^ — > tt'^tt'^X'^Vt, where 
X"^ is the unidentified 4''' track), we require 0.3 GeV 
< M„±„± < Mj to ensure that the charged pions are 
consistent with coming from a t lepton decay. For the 
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Figure 3: The cosine of the angle between the muon and the closest 
identified photon (co.v((?^y)) with all other selection criteria applied for 
about 15% of the BAB/'R data sample. The points represent the data, 
the empty histogram represents the 3 prong r decays, the light shaded 
histogram represents the other t^t" MC, the medium dark histogram 
represent the jU^/i" MC and the dark histogram represents the qq MC. 
The background contamination in these samples is small. 



"opposite-sign" channel (t=^ — > ji'^ji^X'^Vt), we require 
\M„±„T - Mp\ < 100 MeV to ensure that the charged pi- 
ons are consistent with coming from a p meson. This 
produces a loose selection for the "same-sign" channel 
and a tight selection for the "opposite-sign" channel. An 
event can be selected in either or both channels. In the 
case where more than one same-sign or opposite-sign 
pion pairing is possible, the pair with the highest labo- 
ratory frame p, is selected. 

To remove qq backgrounds, events with neutral parti- 
cles with an energy greater than 0.5 GeV that are within 
90 degrees from the muon track are removed. Figure [3] 
shows the cosine of the angle between the muon and the 
photon {cos{9fty)). To suppress radiative di-muon and 
Bhabha backgrounds with conversions, the muon track 
must have a CM momentum, less than 80% and 

greater than 20% of ^/s/2, where ^/s is the beam CM 
energy. To further reduce the non-r backgrounds, the 
polar angle of the system of charged particles, the fi-nn 
system, in the CM frame must satisfy | cos(0ij-„„)\ < 0.8, 
with the net transverse momentum of the fx-Tm system 
being more than 0.3 GeV. 

After the same-sign and opposite-sign events have 
been selected, fourth track candidates are selected, 
which are required to have the appropriate charge to 
come from a r pair event and satisfy the track defini- 
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Figure 4: The 2p^*'/ ^^s of the tag track with all other selection cri- 
teria applied for about 15% of the BABAR data sample. Contamination 
from di-muon and Bhabha events, which peak at 2P/ ■\/s=l.O, are neg- 
ligible. The points represent the data, the empty histogram represents 
the 3 prong r decays, the light shaded histogram represents the other 
r^T'^ MC, the medium dark histogram represent the p^jU" MC and the 
dark histogram represents the qq MC. 
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Figure 5: The track multiplicity in events that have been selected with 
the same-sign or opposite-sign selection presented using the CT and 
GT definitions of the fourth track with all criteria applied for about 
15% of the BABj>R data sample. The points represent the data; the 
contributions from dift'erent backgrounds are shown in the histograms. 
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tions being studied. Figure |5] shows the multiplicity of 
the selected same-sign and opposite-sign events for the 
CT and GT definitions. Once the fourth track candidates 
have been selected, the tracking efficiency is determined 
by using Eq. [T] The difference in the tracking efficiency 
between data and MC is defined using Eq. |2] 



1 



(2) 



305 Similarly, the charge asymmetry of the tracking effi- 

306 ciency is defined using Eq. |3] 

a± = . (3) 

e+ -I- e- 

307 where the efficiency measurements in Eq.|2]and Eq.[3] 

308 also include the detector acceptance. 

309 Monte Carlo studies indicate that the backgrounds 

310 that could potentially bias the determination of the rel- 

311 ative tracking efficiency and the charge asymmetry are: 

312 events with two primary tracks from the e^e" collision 

313 and a photon that converts into an electron pair; qq and 

314 T pair events with six tracks; and t" — > n^K^^Vr where 

315 the decays into a tt^tt^ pair with a vertex that de- 

316 viates significantly from the primary vertex. For the 

317 background events with conversions and /T,, the re- 

318 construction efficiency could differ from that of tracks 

319 originating from the interaction point of the e^e" col- 

320 lision. The largest source of conversions comes from 

321 hadronic t decays with one charged track and 1 or 

322 more neutral particles. This includes — > p'^Vr and 

323 — > h'^nPn^VT (h = TT or K), which have branching 

324 fractions of (25.51 + 0.09)% and (9.51 + 0.1 1)% 10] re- 

325 spectively. The contribution from the r decays with a A'j 

326 is small due to the suppression by the selection cuts and 

327 the branching fractions. The largest background from 

328 events with six tracks originating from the e^e" colli- 

329 sion is from fi^VftVr, — > h*h*li^lT^k^ > 

330 neutrals VT{excludingK^) events which have a branch- 

331 ing fraction of (0.102 + 0.004)%. The contamination 

332 from T pair events with six tracks is <*: 0.1% for both 

333 the same-sign and opposite-sign channels. 
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4.3. Systematic Uncertainties 

The primary systematic uncertainties in measuring 38i 

the tracking efficiency and charge asymmetry arise due 332 

to mis-modeling of background contamination, which 333 

can bias the tracking efficiency due to fake tracks. The 334 

largest background comes from events with two tracks 335 

and a photon that converts in the detector material pro- 386 
ducing an e^e" pair This background includes contri- 
butions from T pair events, radiative di-muon, Bhabha 
events, and two-photon events. 



We estimate the effect of background mis-modeling 
on the efficiency measurement using control samples se- 
lected to be enriched in photon conversion backgrounds. 
The control samples are selected using the standard se- 
lection, minus the vertex requirements, the loose elec- 
tron rejection using PID, and the same-sign and op- 
posite sign invariant mass cuts. Instead of these we 
apply a tight electron PID selection to two oppositely 
charged tracks. The invariant mass of the two oppo- 
sitely charged tracks is required to be less than 0. 1 GeV 
using an electron mass hypothesis. The agreement be- 
tween the data and MC for the selection efficiency of 
this control sample is taken as the uncertainty on the 
modeling of conversions. This is propagated to the A 
and the charge asymmetry measurements using the mea- 
sured rates. Note that this systematic error includes 
both contributions from the mis-modeling of the con- 
versions, and contributions from backgrounds that are 
not included in the MC simulation. 

To assess the impact of potentially different track 
multiplicity from qq backgrounds and the small contri- 
bution from T decays with a K^^, the efficiency differ- 
ence A and the charge asymmetry are calculated with- 
out subtracting these backgrounds. The difference be- 
tween these and the nominal values computed after 
background subtraction is conservatively taken as the 
systematic uncertainty. 

To account for possible differences in the rate of fake 
tracks, a systematic uncertainty based on the difference 
between A and the charge asymmetry calculated with 
(eA) and 

e'xA = ^ (4) 

is included, where A^5 is the number of events where two 
candidate fourth tracks are found. 

As a cross check on the systematic errors, we com- 
pute A and the charge asymmetry, a±, separately in the 
same-sign and opposite-sign channels, and find these to 
be consistent within statistical and background uncer- 
tainties. 

In general, tracks selected in an analysis will not have 
the same kinematic distributions as the tracks in the 
Tau31 study. Therefore, when applying the efficiency 
results of the Tau3 1 study to an analysis, an additional 
systematic uncertainty is needed to account for the ef- 
ficiency dependence on track kinematics. In the Tau31 
analysis we do not estimate the dependence of tracking 
efficiency on track density or track multiplicity. That is 
done in the ISR analysis presented in Section |5] 

We quantify the kinematic variation in A and the 
charge asymmetry by measuring them as a function of 
fourth track polar angle and transverse momentum pt. 
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Figure 6: The cos{9^,h Track) as a function of the cos(6""^'^) for data 
events in Runs 1-6 selected with the same-sign and opposite-sign se- 
lection criteria. The fourth track is identified using the CT definition. 
The dotted lines indicate the boundaries of the coi(6""'') regions se- 
lected for determining the systematic uncertainty on A and the charge 
asymmetry as a function coi(04,;, Track)- 



Figure 7: The Pf,4r/, Track as a function of the p™-" for data events 
in Runs 1-6 selected with the same-sign and opposite-sign selection 
criteria. The fourth track is identified using the CT definition. The 
dotted lines indicate the boundaries of the pj'"" regions selected for 
determining the systematic uncertainty on A and the charge asymme- 
try as a function P,Ath Track ■ 



Because of the three missing neutrinos in the event, 6 400 
and p, of the fourth track cannot be exactly determined. 401 
We therefore construct estimators based on the trajec- «2 
tories of the muon and two identified pions, and use 
them to define different kinematic regions. We define 
the cos{ff) estimator to be 

cos{0""-'-')^cos{e„,„,), (5) 

where the 7Tin2 system is defined as the vector sum of 
the two identified pions. The correlation between the 
coi(0'"''') estimator and the cos(04,/, jmck) is shown in 
Figure |6] 

For p, we define the estimator as 



where yfs is the beam energy, E„j is the energy of the 
ith identified pion and m,r is the mass of a pion. The 
correlation between the p"'""' and P, 4,/, Track is shown in 

404 

405 

The systematic uncertainty on A and the charge asym- 
metry as a function of the estimated cos{d) and p, is 

, 407 

defined as the RMS= Jji" (A/fl± - A,/fl±,,)^ / (« - 1), 408 
where n is the number of regions, A/fl± is the average A 409 
or charge asymmetry as defined previously, and Ai/a-tj 410 



is the A or the charge asymmetry in the ;th region se- 
lected with the estimator The systematic uncertainty 
due to P, and 6 dependence is quantified in Table [1] 



Table 1: Systematic uncertainties for A in P, and 6. The P, and 6 
uncertainty in the T(2s) and T(3s) runs are sensitive to the limited 
statistics. 



Data Period f , Uncertainty 


Uncertainty 


GTL Track Definition 




Runs 1-6 0.20% 


0.11% 


Run 7 - r(2s) 0.30% 


0.38% 


Run 7 - r(3s) 0.80% 


0.42% 


CT Track Definition 




Runs 1-6 0.10% 


0.06% 


Run 7 - r(2s) 0.65% 


0.28% 


Run 7 - r(3s) 0.79% 


0.32% 



4.4. Tau31 Results vs. Run Period 

Figure[8]shows the run-by-run tracking efficiency for 
the two track definitions studied in this analysis: GT 
and CT. The tracking efficiencies in the data and MC are 
found to be consistent with each other This can be seen 
in Figure |9] which presents A. The charge asymmetry 
can be seen in Figure [TOl The plots suggest that there is 
no significant charge bias in the tracking efficiency. 
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Figure 8: The tracking efficiency as a function of run number for the 
GT definition (top) and the CT definition (bottom). The data are repre- 
sented by solid markers and the MC by open markers. The error bars 
represent the total uncertainty where the statistical and background 
systematic uncertainties have been added in quadrature. 



The stability of the agreement between data and MC 
over the 7 run periods (as shown in Figures |9] and [TOl l 
demonstrates that the detector simulation, which is up- 
dated regularly, accurately models the tracking perfor- 
mance of the detector as a function of time. Because 
there is no significant time variation observed between 
Runs 1 and 6 in A and in the charge asymmetry of the 
tracking efficiencies, an average of Runs 1-6 for A and 
the charge asymmetry of the tracking efficiencies is cal- 
culated. These averages are used to calculate the sys- 
tematic uncertainty due to tracking efficiency. The sys- 
tematic uncertainty per track for a given track definition 
is 



^Tracking 



C &CT/GT 

1-Ar 



(7) 



^CT/CT 

where ctact/ct is the total uncertainty on A for the given 
track definition. These results are the primary source of 
systematic uncertainty in track reconstruction efficiency 
in BaBar. 



428 5. Tracking efficiency using the ISR channel 

429 71^ 7l~ 71^ 7l~ y ISR 

430 A complementary approach to the Tau31 method is 

431 to study the tracking efficiency using processes such as 

432 e^e" — > TT^n^yisR and e^e^ n'^n^n^n^jisR, where a 

433 high energetic photon jisrK emitted from an initial lep- 

434 ton. This final state provides a clean event sample, cov- 

435 ering a wide range of momenta and polar angles of the 

436 tracks. In this section, we describe one such measure- 

437 ment involving four pions in the final state along with 
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Figure 9: The data-MC difference in the tracking efficiency as a func- 
tion of run number for the GT definition and the CT definition. The 
error bars represent the total uncertainty where the statistical and sys- 
tematic uncertainties have been added in quadrature. 



e 2 

1 


-1 

-2r 
2: 
1 

Or 
-1r 
-2 



■ GTData 
□ GTMC 



• GTData 
o GTMC 
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Run 1 Run 2 Run 3 Run 4 Run 5 Run 6 Run 7 

Figure 10: The charge asymmetry of the tracking efficiency as a func- 
tion of mn number for the GT definition (top) and the CT definition 
(bottom). The data ai'e represented by the solid mai'kers and the MC 
by the open markers. The error bars represent the total uncertainty 
where the statistical and systematic uncertainties have been added in 
quadrature. 
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438 the ISR photon. The Tau3 1 method has a higher statist!- 484 

439 cal accuracy, allowing the explicit study of time depen- 485 

440 dent effects. By contrast, since no neutrinos are present 486 

441 in the final state, the ISR events allow a more precise 487 

442 estimate of the missing track parameters than the Tau3 1 488 

443 method. In addition, the track density for ISR events is 489 

444 higher compared to the events in the Tau3 1 study, cor- 490 

445 responding to different B4R4R physics channels. The 491 

446 high track density in combination with the precise track 492 

447 parameter prediction allows studying the track overlap 493 

448 effects in tracking efficiency. 494 

449 To study tracking efficiency with ISR, we use two 495 

450 event samples: one in which all 4 charged particles 

451 are reconstructed (4-track), and one in which only 3 

452 charged particles are found (3-track). Using energy and 

453 momentum conservation in a kinematic fit, we can accu- 

454 rately predict the direction and momentum of the miss- i98 

455 ing track in the 3-track sample. By calculating the ra- 199 

456 tio of the number of lost tracks Niosi tracks to the number 5™ 

457 of measured tracks Ndeieaed tracks, we obtain the tracking 

458 inefficiency, defined in equation (O, and the track- ™2 

459 ing efficiency, e, according to equation (|9]). Both can be 

460 measured as a function of the kinematic properties of 

461 the missing track. 505 



tracks 



^detected tracks tracks 



(8) 



(9) 510 



5. 7. ISR Event Selection 



463 For the ISR efficiency measurement we require that 513 

464 the tracks have a polar angle inside the detector accep- 514 

465 tance (-0.82 < cos 9ch < 0.92), and that the transverse 515 

466 distance of closest approach of the track to the event 516 

467 vertex (or nominal interaction point if no primary event 517 

468 vertex is found) be smaller than 1.5 cm, and be within 518 

469 2.5 cm in the beam direction. Tracks with less than 519 

470 100 MeV/c transverse momentum are rejected. The ISR 520 

471 photon is restricted to the polar angular range inside the 521 

472 EMC acceptance (0.5 rad < 0y^^^ < 2.4 rad), and a mini- 522 

473 mum photon energy of Eisr > 3 GeV is required. Either 523 

474 3 or 4 selected tracks are required in the event. 524 

475 In order to suppress radiative Bhabha events, we re- 525 

476 ject events where the two most energetic tracks pass a 526 

477 loose electron PID selection. This also removes most yy 527 

478 events with an additional high energetic photon {Ey^cm > 

479 4 GeV) in opposite direction to the ISR photon can- 529 

480 didate. We require the minimum angle between the 530 

481 charged tracks and the ISR photon to be larger than 531 

482 1.0 rad, which rejects a large fraction of e^e^ — > qq, 532 

483 (q - u, ,s) and e^e^ — > r^r" event backgrounds. Events 533 



with one or two tracks with PID consistent with being 
a K'^ in the 3-track or the 4-track sample are rejected, 
respectively. Finally, we require the 4-n invariant mass 
to be in the range of l.lGeY/c^ < M4„ < lAGeV/c^, 
where we expect a high signal to noise ratio. 

Backgrounds from e^e^ qq (q = u, ,s) are sim- 
ulated with JETSET fl^l, while e+e^ -» t+t^ back- 
grounds are simulated using KORALB 111 If . The ISR- 
channels are simulated with the AFKQED 112 [genera- 
tor, based on an early version of PHOKHARA 01311 . The 
MC samples are normalized according to the luminosity 
observed in the data. 



5.2. ISR Kinematic Fit 

Selected events are subjected to a kinematic fit as- 
suming the TT^TT^n^n^yjsR signal hypothesis;if4^, as well 
as the K^K^TT^TT^yisR background hypothesis X2K2n- 
The fit in the 4-track (3-track) sample uses the four 
(three) tracks, the ISR photon and the kinematic infor- 
mation of incoming electron and positron. Energy and 
momentum conservation leads to four (one) constraints, 
or a 4C-fit (IC-fit), respectively. 

The resulting x~ distributions are shown in Fig. fTTl 
The distributions are broader than expected, partly 
due to detector resolution effects, but mostly because 
additional ISR photons are not included into the kine- 
matic fit hypothesis. In Fig. [TTla) the 4-track sample 
shows a good agreement between the data and MC in 
the presence of negligible background. 

In Fig. fTTT b) the corresponding distributions are 
shown for the 3-track sample. Here, we also require 
the predicted polar angle for the missing track be in the 
detector acceptance (-0.82 < cos 9ei, < 0.92) region. 
The relative amount of background is much larger in 
this sample, since the kinematic closure that suppresses 
a lot of background in the 4-track sample is weaker with 
only one constraint. The visible difference between the 
number of events after background subtraction (red) and 
signal MC (black) suggests that more n tracks are lost 
in data than are described by MC. 

Fig. fTTT b) also shows a difference in the shape of the 
X^ distributions between the data and MC. The plateau 
in the background MC at large suggests that all back- 
grounds are not subtracted from data. Therefore we 
perform an additional background subtraction based on 
data sidebands. The idea of the sideband subtraction is 
illustrated in Fig. [121 which plots the 3-track distri- 
bution for a subset of the BaBar data. We define a sig- 
nal region enriched in signal events, which contains A^i 
events. The control region, which has substantial back- 
ground contributions, contains N2 events. Let A^i , (Nib) 
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Figure 11: (a): x^^^ distribution for 4-track sample (4C). Data with 
subtracted background (red points), signal MC (black histogram) and 
the sum of background MC channels (blue histogram), (b): Corre- 
sponding distributions for the 3-track sample (IC). The signal and 
control regions are indicated with vertical lines, with the region in the 
extreme left being the signal region and the area in the middle being 
the control region. 



be the number of signal (background) events in the sig- 
nal region, and N2s {N2b) the corresponding numbers for 
the control region. Assuming one knows the ratios, 

N2h 



a = and b — (10) 



the number of signal events can then be calculated ac- 
cording to the following equation: 



b-Nx-N2 



(11) 



We define signal and control regions in the 4-track ses 
sample as x\„ 4c < 30 and 30 < xl,^ 4c respec- see 

tively. The corresponding regions in the 3-track sample se? 
are chosen so that the ratio of events in the signal to con- sea 
trol region is the same as in the 4-track sample, resulting seg 



< 3 and 3 < ic ^ ^ respectively. The 



tio a is determined using signal MC. In order to obtain 
b, we assume any difference in tracking inefficiency be- 
tween data and MC does not depend on;^'^^. Therefore 
we performed a fit of the difference between data and 
MC using a linear Probability Density Function (PDF), 
allowing a scale-factor for MC. 
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Figure 12: Fit result for sideband parameter h using Run 5 fitting 
signal MC (blue histogram) and a linear background (blue line) to 
data (black points). Also indicated are the number of signal N is 12s 
and background Ni),/2b events in the signal and control region. 

The result of the fit is shown in Fig. [12] Small dis- 
crepancies at low x^ indicate that there is still some 
background present. The remaining difference in the;^f^ 
distribution is consistent within the uncertainty of the 
cross section of the peaking background contributions 
that have been subtracted. After subtracting the addi- 
tional background using equation [TT] the inefficiency 
difference between data and MC A77 - rjdaia - Hmc is 
determined to be 



At] = (0.75 ± 0.05,,,,, + 0.34,,,,,)%. 



(12) 



The systematic uncertainty on At] is dominated by 
the uncertainty of the cross section of the subtracted in- 
dividual background contributions in the 3-track sam- 
ple. Most of these cross sections have been measured 
in previous BaBar analyses 

idEHInllil. The 

normalization of the additional contributions of contin- 
uum and e^e" — > tt backgrounds have been verified 
with specific kinematic distributions. Note that this re- 
sult is not directly comparable to the Tau31 efficiency 
result, as that was calculated using an isolation require- 
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570 ment between the tracks. The effect of track overlaps is 

571 discussed in the next section. 

572 5.3. ISR Efficiency Kinematic Dependence 

573 In Fig. [13] the dependence of At/ on the polar angle 

574 (a) and the transverse momentum p, (b) of the missing 

575 track is presented. The dependence on p, is flat within 

576 the uncertainties of 0.4%. A slight dependence on the 

577 polar angle is visible with almost no difference between 

578 data and MC in the forward region at small polar angles 

579 and a difference of approximately 1 % in the central and 

580 backward region. Due to the beam energy asymmetry 

581 at BaB/>R, high energy photons are preferably emitted 
682 in the forward direction at small polar angles. In ISR 

583 events, the hadronic system is emitted back-to-back to 

584 the ISR photon. The energy of the photon is correlated 

585 with the opening angle of the cone of the hadronic sys- 

586 tem. This correlation leads to an increasing track over- 

587 lap probability in the backward region of the detector, 

588 which is not perfectly modeled by MC as shown in the 

589 following. 

590 One source of tracking inefficiency is when two 

591 tracks overlap in the detector, causing sensor signals 

592 from one or both to be lost or distorted, and creating 

593 hit patterns that can be hard for the track finding al- 

594 gorithms to distinguish. BaBar tracking inefficiency 

595 is most affected by overlaps in azimuth, as the DCH 

596 largely projects out track polar angle. Due to magnetic 

597 bending, tracks with the same charge are more likely 

598 to overlap in azimuth than tracks with opposite charge. 

599 Furthermore, the overlap between tracks with opposite 

600 charge depends in an asymmetric way on the azimuthal 

601 angle between them. These effects are shown schemat- 

602 ically in Fig. [141 To study the dependence of tracking 

603 efficiency on overlap, we define variables sensitive to 

604 the charge-dependent two-track azimuthal separation: 

605 The effect of track loss due to overlapping tracks with 

606 different charge (DC) is visible in the distribution of 

607 the charge oriented azimuthal angle difference between 

608 the lost track and the reconstructed track with different 

609 charge Acpoc - (pi^^) " <P(^ )- Since in our study there 

610 are always two pions with the different charge of the 

611 missing pion, two angles are obtained for each event. 

612 In Fig. [TsT a) the Acpoc distribution is plotted for data 

613 and MC. The asymmetric distribution shows that the DC 

614 tracking inefficiency peaks at small positive values. 

615 The number of tracks lost due to DC track over- 

616 lap is estimated by subtracting the negative half of this 

617 distribution from the positive, as illustrated for MC in 

618 Fig. [IST b). The inefficiency is corrected as indicated 

619 in equation [13] leading to a correction for overlapping 

620 tracks with different charge DC of 0.41 %. 




-0.005 - 

-o.oi l" 

0.5 1 1.5 2 2.5 3 

p,[GeV/c] 

Figure 13: Relative data-MC dilference of tracking inefficiency vs. 
the polar angle of the track 8 (a) and vs. the transverse momentum 
Pt (b). Red lines indicate the detection region used to determine the 
average inefficiency. 



, ^ lost tracks ^ overlapping tracks ,-,^s 

rj - (13) 

A^ tracks 

621 We describe the same charge (SC) track overlap in- 

622 efficiency in terms of A0sc = \<Pi''^^) — <Pi^^)\^ illus- 

623 trated in Fig. [14] (b): the angle between the lost track 

624 and the reconstructed track with the same charge. For 

625 data in Fig. [TsT c) the angle between lost track and re- 

626 constructed track with the same charge in the 3 -track 

627 sample is plotted in red. The blue histogram shows the 

628 same distribution for the two detected tracks. The dis- 

629 tribution with one lost track is the superposition of the 

630 distribution due to detection inefficiency and a peaking 

631 distribution at small A4>sc due to track overlap losses. 

632 The distribution due to usual detection inefficiency has 

633 the same Atpsc dependence as the distribution of the two 

634 measured tracks. The number of tracks lost due to track 

635 overlap can be estimated by scaling down the distribu- 
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Figure 14: (a): Charge oriented azimuthal angle between missing pion 
and detected pions with dilferent charge A0oc = <f'i^*) ~ <f>{^~) (2 en- 
tries per event); (b) absolute value of azimuthal angle between missing 
pion and detected pion with same charge A(psc = lipi^*) ~ 



tion of the measured tracks until the tails of the distri- 
bution match with the distribution including one miss- 
ing pion. The difference at small Aipsc is a good esti- 
mate for the number of tracks lost due to track overlap. 
The corresponding distributions for MC are displayed 
in Fig. [TSl d). The effect of SC tracks overlap is well 
modeled in MC. 



5.4. ISR Efficiency Summary 

To summarize, the difference in tracking inefficiency 
per track including track overlap is determined from 
ISR events to be: 



At? = (0.75 ± 0.05 + 0.34,,,,)% 



(14) 



Because of the track isolation requirement applied in the 
Tau3 1 selection, the different track multiplicity, and the 
different event topology, the ISR study includes a sig- 
nificantly higher track overlap probability and thus the 
value in equation[T4]is not directly comparable with the 
Tau3 1 result discussed in Section|4] To make a compari- 
son, we quantify the effect due to track overlap by study- 
ing the distributions of the azimuthal angular difference 
between same charged tracks and oppositely charged 
tracks. Taking this effect into account, we measure an 
efficiency difference between data and simulation of: 
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Figure 15: (a): Angle between missing pion and detected pions with 
different charge for the data (red points) and signal simulation normal- 
ized to the luminosity (black histogram). Two entries per event, (b): 
Illustration of the cleaning procedure. Tracks lost due to DC overlap 
(red) and due to other effects (black), (c): Angle between missing 
pion and detected pions with same charge (red) and the angle between 
two detected pions (blue) normalized to the same number of events in 
the region 0.3 rad < A<I) < 0.8 rad. (d): Same as (c), but for signal 
MC. 



At?' = (0.38 + 0.05 + 0.39,,.,,)% 



(15) 



658 This result is consistent with the Tau31 efficiency 

659 difference within the uncertainties. Depending on the 

660 event multiplicity and kinematics, BaBar analyses may 

661 need the inefficiency with or without track overlap ef- 

662 fects. 



6. Tracking Charge Asymmetry 

664 Since a main objective of the BaBar experiment is to 

665 measure CP violation, it is vital to understand and mea- 

666 sure any possible charge asymmetry in the track recon- 

667 struction. For instance, a promising mode for searching 

668 for CP violation in charm decays is D'^ — > K^K^ii^. 

669 An asymmetry in the reconstruction efficiency for the 

670 n"^ would bias the CP result. Because the signal in these 

671 decays has a statistical uncertainty of ~ 0.25%, a com- 

672 parable control of the tracking efficiency asymmetry is 

673 needed. 

We define the charged pion tracking asymmetry as 



a{pLah) = 



e{p„+) + e{p„-) 



(16) 
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Figure 16: Tracking asymmetry from MC as a function of charged 
pion momentum in the laboratory frame. The inset plot shows the 
asymmetry up to 2 GeV. 



where indicates that momenta are in the lab frame, 
and p„t (p„-) refers to the momentum of the positively 
(negatively) charged pion. 

We illustrate our expectations in this regard using 
MC. Figure[T6]shows the pion tracking efficiency asym- 
metry derived from MC using generator information 
for pion tracks in D'^ K^K^ji^ decays. The av- 
erage asymmetry for MC in this mode is found to be 
a(pLab) = (-6 + 23) X 10"^, consistent with zero within 
the uncertainties, and without any significant momen- 
tum dependence. 

Two different methods are used to determine the pion 
track efficiency asymmetry directly from data. The 
more precise technique relies on Tau3 1 events. We work 
directly in the observed variables and use the ratios of 
the numbers of two-hadron decays to three-hadron de- 
cays to determine the pion inefficiency. Instead of fitting 
distributions of 2- and 3-hadron decays, we recognize 
that the (fewer) 2-hadron events that arise from track- 
ing inefficiency can be easily modeled directly from the 
3-hadron events. In practice this is done by weighting 
every 3-hadron event by the ratio (1 - e)/e, where e 
is the track efficiency of the observed third track. For 
both 3-hadron as well as 2-hadron events we select only 
events from the p-decay channels — > p'Vj^v^, ac- 
cording to the selection criteria described in section |4] 
since the inclusive — » n^n^h^Vj has more signifi- 
cant backgrounds, specifically with contamination from 
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Figure 17: Distribution of the observed 2-hadron events for the full 
B4R4R data sample. Bin boundaries were chosen to obtain the same 
number of events (215) in each bin. 



electrons. The total number of 2-prong (3-prong) events 
in the sample is 86,092 (1,365,900). The distribution of 
events in the observed variables, pf„„sj and cos(0)„,;.55, is 
shown in Figure [17] The observed variables are deter- 
mined from the 2-prong momenta: 



p(nn) = p(n*) + p(n ) 



such that 



and 



pinn) 



(17) 
(18) 
(19) 



In order to fit these event distributions, one must 
also account for backgrounds. The 2-hadron events of 
interest include approximately 7% background events. 
Chief among these are events from photon (5%) and 
tt" (1%) conversions in 1-hadron decays of the tau, 
where the 1 -prong track from the tau is combined with 
a track from the photon or tt" and identified as a 2-prong 
event. Inelastic nuclear interactions due to tracks pass- 
ing through detector material and other backgrounds 
are small in comparison. The backgrounds are split 
into "photon" and "other" components and the overall 
normalization of each distribution is a parameter in a 
binned fit. Another large background contribution to 
2-hadron events (whose normalization is a parameter) is 
acceptance loss events due to the third track being lost in 
the direction of the beam. PDFs are obtained from MC 
as normalized histograms in the observed variables of 
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the various backgrounds; events in these have been re- 
weighted to account for inadequacies in the MC 3-body 
Dalitz distributions by matching the 3-body mass distri- 
bution as well as both the 2-body mass distributions to 
those in data. The tracking efficiency asymmetry fit is a 
binned ;^f^-fit with binning as shown in Fig.fTTl 

The significant parameters in the fit describe the 
tracking efficiency and the asymmetry as a function of 
the lab momentum. The tracking efficiency is parame- 
terized with the following phenomenological formula: 



P Lab- Pi) 
^0 



e(pLah) = 1 - Aoe ^0 - Bqc n (20) 

where the parameters are Aq, Bq, pq, pi, tq, and t\, in 
addition to parameters which measure the asymmetry in 
bins of lab momentum. Finally, it should be mentioned 
that we account for differences in the 3-hadron distri- 
butions of OTj2 versus (1,2,3 denote the particles 
in the 3 -prong tau decay) in data and MC by weighting 
3-hadron events according to the data/MC , m^^ dis- 
tribution ratio. The fit to our data is good as evidenced 
by a;i'VNDF = 792/780, i.e., a 37% probabihty. 

Results from this procedure are shown in Figures [18] 
and [19] We find the average charged pion tracking effi- 
ciency asymmetry to be (^(pio/,) = (0.10±0.26)%,inour 
momentum range of approximately 0-4 GeV/c, consis- 
tent with zero. To account for systematic errors we re-fit 
the data with the following variations. We force the ac- 
ceptance loss and background descriptions in the fit in- 
dividually to be charge-independent, and we reduce the 
number of background components by combining some 
PDFs. We find the total systematic error to be 0. 10%. 
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Figure 18: The tracking efficiency determined by the Tau31 method 
as a function of charged pion momentum in the laboratory frame. The 
red envelope around the efficiency curve indicates Itr statistical error 
bands. 

Another technique we use to measure the charged 
track efficiency asymmetry utilizes isotropy of spinless- 
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Figure 19: The tracking asymmetry determined by the TauSl method 
as a function of charged pion momentum in the laboratory frame. The 
average asymmetry over momenta 0-2 GeV/c determined from D" 
decays is also shown here for comparison. 



two-body decays. In this method we study the D — » 

—0 

n n and D — > n n decays. We require that these de- 
cays not be from B-meson decays (as these have larger 
backgrounds) and that they be tagged as being from D*^ 
decays to improve signal purity. Also, in both cases we 
require that at least one pion have momentum greater 
than 2 GeV/c and assume that the tracking efficiency 
charge asymmetry is zero for this pion. Therefore, any 
asymmetry in yields is the result of tracking asymmetry 
of the lower momentum pion which is reported below. 

n — 

High purity samples of D and D decays are ob- 
tained using slow pions associated with the decay of D*^ 
to tag the flavor of the meson. A detailed description 
of the event selection is described in the publication of 
£)0 _ ffi mixing using the ratio of lifetimes for the decay 
Particle identification is not applied 
to the selection of pion tracks, rather we choose to re- 
move reflections from the K^n^ decays of the using 
a cut on the reflected mass and we account for the re- 
maining contamination from the tails by studying their 
Ti^n'' mass distributions and including a term with such 
a shape in our 1-D binned n^n'' mass fit. Yields of D" 
decays where the higher momentum track is either the 
71^ or the ;7r" are separately determined and are used to 
determine the asymmetry. A similar study is carried out 
using D decays, and the combined charge asymmetry 
of the efficiency, averaged over pion momenta from to 
2 GeV/c is found to be a{pLah) = (-0.12 + 0.50)%, con- 
sistent with zero and the Tau3 1 method result, but not 
as precise as the Tau3 1 method. 
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7. Low pt tracking efficiency measurement 

The T pair sample provides an estimate of tracking ef- 
ficiency for charged tracks with pj- > 180 MeV/c only. 
However, the detection of low pj tracks {pj < 180 
MeV/c) is important for tagged D° analyses. tagging 
is performed through the D*^ — > D" n'^ decay, where the 
soft pion {nl.) is emitted with an energy just over its rest 
mass in the D*^ frame, and so typically has very low 
Pt in the lab frame. D" tagging is used in CP violation, 
mixing, and many other precision analyses, therefore a 
good understanding of the low pj tracking efficiency is 
required. 

The low Pj reconstruction efficiency analysis is 
based on a previous analysis by the CLEO collaboration 
ll20ll . CLEO demonstrated that the relative slow pion ef- 
ficiency can be measured as a function of momentum 
using helicity distributions. The slow pion helicity an- 
gle 0* is defined as the angle between the slow pion mo- 
mentum in the D* rest frame and the D* momentum in 
the laboratory frame. This is illustrated in Fig.l20l 




just on the cosfl*, but also on p*{D*), we perform an 
angular efficiency analysis in bins of p*{D*). We then 
fit these cos 0* distrubions to a function defined as the 
convolution of Eq.|2T|and the efficiency function: 



e{p) 



1 



S(p-po)+l ' 



if P > Pq 
ifp < p(). 



(23) 



The goal of this analysis is to compare data and MC 
efficiencies to get a systematic error from the relative 
difference between them: 



/ edala(p)dp - J £Mc(p)dp 
J edata{p)dp 



(24) 



The limitations of this method are the effects that may 
be not correctly described in the MC, such as final state 
interactions or radiative losses. 
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Figure 20: Definition of slow pion helicity angle ( 
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When a vector meson decays to a final state made 
of two pseudo-scalar mesons, the distribution of the he- 
licity angle is expected to be symmetrical and can be 
described as ll2lll22ll 



dN 



d cos 0* 



oc (1 H- acos 0*), 1 < a < -Hoo, 



(21) 



Furthermore, the cosine of the helicity angle is related 
to the slow pion momentum by: 



Pn.. = y{pl^ COS 0* -PEI), 



(22) 



where fi and y are the D* boost parameters. Since p* and 
E* are known once the D* momentum is known, Eq. |22] 
maps any asymmetry observed in Eq. |2T| to a relative 
reconstruction inefficiency in a specific part of the slow 
pion momentum spectrum. 

We measure the cos 0* distribution in 8 bins of p*{D*) 
spectrum as shown in Fig. |2T| Since p„ depends not 



Figure 21: Distribution of p*(D*) in the data sample. On top of the 
figure the blue lines show the lower and upper limit of each bin indi- 
cated by the red number. 

The analysis is done using 470 fb"' of data recorded 
by the BABm detector and about 4.2 x 10^ generic MC 
events. The decay chain e+e" — > D*^ X, D*^ 
TTj , D° — » K JT'^ 12311 is reconstructed in both data and 
MC, requiring particle identification for the kaon and 
the two vertices to be successfully reconstructed. The 
D° — > mode is chosen to provide a clean sample 

of D*^ D° TT^ decays with a high branching fraction. 
A control sample is reconstructed the same way by not 
requiring the kaon identification. This sample is used 
for background subtraction. The p*{D*) spectrum has 
been compared between data and MC. Differences are 
corrected for by weighting the MC sample which is then 
normalized to data. 

As shown in Fig.|22] four categories of events can be 
recognized after the reconstruction: 
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1 . signal: real D" and n'!, from D*^ decay. 

2. Missed ;r^: real D'^ — > K decay that may or 
may not have come from a combined to a tt'^ 
from combinatoric. 

3. Missed Z)": a mis-reconstructed D" with a real tt^ 
from This is mostly D" ^ /T" D" ^ /T" 



7r+ tt", 



TT^ TT or cases where the kaon and 



pion assignments have been swapped. 
4. Combinatoric background: neither D'^ or tt^ are 
correctly reconstructed from a D*+ decay. 

The amount of combinatoric and real D", fake n'^ 
background in the signal region is estimated using the 
re-normalized distribution of the control sample 
sidebands in the Am - m{K^n^n^) - m{K'n^) signal 
region. The scale factor needed for going from Am side- 
band to the Am signal region is taken from the control 
sample itself. This background subtraction procedure 
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Figure 22: m(K~ Tt*) vs. Am scatter plot of the data sample. Signal re- 
gion is identified by the red box and the blue Unes show the sidebands. 
The numbers identify the events category. 



830 has been carried out for the cos 6* distribution for each 

831 bin of p*{D*) using the following steps: 



856 
857 
858 
859 



862 
863 
864 
865 
866 
867 



1 . consider the no PID sample in the m{K^n^) side- 
band regions; divide the number of events in the 
Am signal region by the number of events in the 
Am sideband to get the scale factor to go from side- 
band to signal in Am; 

2. in the good PID sample scale the m{K^n^) spec- 
trum in the Am sideband region using the factor 
obtained in the previous step; then integrate the 
resulting n\{K~n*) spectrum to get the factor to 
rescale background; 

3. use the factor measured in step 2 to rescale 
the interesting distribution (cos (f ) obtained from 



the events of the no PID sample in Am signal, 
m(K^n^) sideband region (category 4); 
4. subtract the distribution obtained in step 3 from the 
same distribution retrieved from good PID sample 
in signal region. 

This procedure has been carried out for the cos 8* dis- 
tribution for each bin of p*(D*). All the histograms have 
been then fit to the convolution of Eqs.|2T|and|23]to de- 
termine the parameters of the efficiency function po and 
S. Measuring pQ and 6 for data and MC, we can evaluate 
the systematic error using Eq. |24] The fit makes use of 
a global defined as 



X 



=2 

l.k 



(D,k-Sitf 



(25) 



where k is the index referring to one of the 8 p*(D*) 
regions, I refers to one of the 16 bins of the cosO* his- 
togram in that region. D^, o-On. and S ik are the number 
of events observed in the bin, its error and the number 
of events expect by the fit model, respectively. The ex- 
pression of the fit model is 



where ; indicates the bin of the cos 0* distribution in the 
j^ih region, and j is one of the 10 bins of the de- 

tailed distribution within the range of momentum con- 
sidered in the k''' p*{D*) region. 

The number of floating parameters in the fit are 18: 
8 normalization factors A^,-, 8 a, (one for each bin) from 
Eq.im and 6 and po from Eq.[23] The fit has been made 
both to data and MC, giving the results shown in Fig.l23] 
and Tab.|2] 

Finally, the efficiency functions are compared in 
Fig- ESI Please note that these distributions include ac- 
ceptance. The method shown herein does not allow to 
disentangle the acceptance from the soft pion efficiency. 
The systematic uncertainty estimated from Eq. |24] for 
the low pt tracks is cr^y.,, - 1.54%. 

8. reconstruction efficiency measurement 

A significant number of analyses in BaBar involve the 
reconstruction of the decay K'-^ — > tt'^tt^, where the two 
charged pions belong to the list CT of all reconstructed 
tracks in the event. The track reconstruction efficiency 
for charged tracks originating within 15 mm in XY from 
the beam spot is studied by the other methods presented 
earlier in the paper However, most of the K'^ 's decay 
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Table 2: Fit results on data and MC. 
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Figure 23: Fit of the model distribution to data (top) and MC (bottom). 
In both the plots, the measured data/MC event ratios are represented 
by the dots, while the fit results are shown using a line histogram. 
The distributions of cos 6* in the different ranges of p'( D' ) are shown 
using different colors, as outlined in the legend (the p*(D') values are 
measured in GeV/c). The fit to data returns ^^/nj^y = 1.34; the one 
to MC gi\es x^/njof = 0.71. 
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Figure 24: Comparison of the fit results on a in the 8 bins of p*(D*) 
for data (black) and Monte Carlo (red). The difference observed in the 
4'*, 5'* and 8''' bins is due to the slightly different helicity distribution 
for data and MC in these p'{D') ranges. 



879 outside this 15 mm radius, making it necessary to under- 

880 stand the K'-^ daughter reconstruction efficiency in data 

881 and MC. 

882 The reconstruction efficiency of the K'^^ daughters de- 

883 pends on the K'^ transverse momentum, pr, polar angle, 

884 0LAB and transverse (XY) flight distance, dxr, which is 

885 computed as the distance between the primary vertex of 

886 the event and the refitted A'" decay vertex. 

887 The general strategy is to subdivide the data and MC 

888 events into a large number of samples by choosing an 

889 appropriate binning in these variables, determine the 

890 number of K'^ 's in each bin in data and MC and, for each 

891 of ffie momentum and polar angle ranges, normalize the 

892 ratio of its value in the first bin in dxr, where all tracking 

893 effects are understood to 1 .000 by definition, with no as- 

894 sociated error other than the systematic uncertainty per 



895 track, as discussed in Section|4] Bin sizes are optimized 

896 to ensure a sufficient number of events in each bin, with 

897 4 bins in /jt- (0.0 - 0.4 - 1 .0 - 2.0 - 4.0 GeV/c), 8 bins in 

898 6lab (7.0° - 25.6° - 44.25° - 62.88° - 81 .5° - 100.13° - 

899 118.75° - 137.38° - 156.0°) and 10 bins in dxr (0.0 - 

900 0.3 - 1.3 - 2.78 - 3.2 - 4.0 - 5.4 - 9.1 - 11.4 - 23.6 - 40.0 

901 cm). The binning in dxr roughly reflects the structure 

902 of the BaBar detector and the bins are numbered from 

903 to 9, bin 1 being the normalization bin. The normaUzed 

904 ratio of the data and MC as a function of different bins 

905 is provided as a correction factor for the K'^ daughter 

906 reconstruction efficiency. In order to reduce the uncer- 

907 tainty from the imperfect simulation of the random track 

908 background or the potential differences between the K'^ 

909 quality cut efficiencies in data and MC, we remove the 
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Figure 25: Soft pion reconstruction efficiency functions obtained from 
the fit to data (red) and MC (blue). Both the efficiency functions are 
shown together with the functions obtained by varying the central val- 
ues of the fit parameters po and 5 by Icr. The curve obtained using the 
central values is drawn in black. On top of the curves, the distribution 
of the relative difference between data and MC is shown. 



910 immediate vicinity of the event's primary vertex which 

911 is 3 mm in XY from the first (normalization) bin. 

912 Events of interest are selected by looking for the 

913 B — » h^h'K'-^ (with h - n,K) decays in the data and 

914 MC samples. The MC sample includes events from 

915 generic B decays, light quark events (u,d,s,c) and t^t" 

916 decays. The K^^ is reconstructed from two oppositely 

917 charged tracks, the invariant mass of which is required 

918 to be within 25 MeV/c^ of the PDG value of the 

919 mass (m^o = 497.614 + 0.024 MeV/c^) The two op- 

920 positely charged tracks must originate from a common 

921 vertex and the fit is required not to fail. The event is re- 

922 quired to have at least five GT tracks, two of which are 

923 oppositely charged GT tracks that when combined with 

924 the candidate to form an object with mes > 5.19 

925 GeV/c^ and lA^I < 0.3 GeV. About 93% of these events 

926 come from the light quark (udsc) continuum; the con- 

927 tribution from r^r" production is about 2.5% and only 

928 about 3.5% of candidates arise from B decays, most of 

929 which are random track combinations. Figures [26] and 

930 122] show the data and MC comparison of the A'^' mass, 

931 pt, Slab and dxY distributions for the reconstructed A'^' 

932 candidates. 

933 To determine the number of K'^'s in each of the bins 

934 in data and MC, the A"" mass distributions in each of 

935 the bins are fitted with a sum of a double Gaussian and 

936 a constant background. The constant background, de- 

937 termined from the sideband regions in the mass dis- 

938 ti-ibution, [0.476, 0.485] U [0.511, 0.520] GeV/c", for 

939 each bin, is then subtracted to determine the number of 
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Figure 26: mass (top) and transverse momentum (bottom) for data 
and MC. MC is normalized to the data luminosity. 



in each bin. The limited statistics in a large fraction 
of bins makes it impractical to allow the slope of the 
background to float. The binned maximum log likeli- 
hood method is used since the default ;(f^-minimization 
method is less appropriate in this study as it systemati- 
cally, in a statistics-dependent way, underestimates the 
number of events in each bin. 

We define the values of the normalized ratios Rijit and 
the uncertainties 0-^^^^, where the indices /, j and k stand 
for the number of the pr, Olab and dxr bins respectively, 
to be 



Rijk^(Nijk/Mijk)/{Niji/Miji) 



(27) 
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Figure 27: K'^ polar angle (top) and transverse flight length (bottom) 
for data and MC. MC is normalized to the data luminosity. 



(28) 



where Nijk and a^, 
uncertainties in data and Mijk and cr^.^^ 



, are the numbers of /T^'s and their 992 
are the numbers 993 
of 's and their uncertainties in MC. We also take into 994 
account the difference in the mass resolutions in data 995 
and MC by performing numerical integration to deter- 996 
mine the efficiencies of the \m„*,j — to^| < 12 MeV/c^ 997 
cut in data and MC. The efficiency corrected normalized 998 
ratio, Rijk, also called the correction factor, is computed 999 
for each bin and also as a function of several different 1000 
sets of K^^ quality cuts (cuts on the K^^ mass, 3-D or XY 1001 
flight length or its significance, the 3-D or XY angle, 1002 



a between the A'^ momentum and the line connecting 
the A"" decay vertex and the primary vertex of the event, 
and, in one instance, a cut on the AT" — > n^ji^ vertexing 
quality), which also enables us to study the systematic 
uncertainties associated with a correction factor 

To get the overall correction in the K^^ daughter re- 
construction efficiency in an analysis, we apply the cor- 
rection factors to the signal MC. If Hjjk is the num- 
ber of events in the signal MC sample that falls within 
the bin (ijk), the relative weight of a K^^ reconstruc- 
tion efficiency correction table element Rijk is Hijk / 
Htot, where - EHijk and its statistical uncertainty is 
'\J{Hijif)/Htot- The central value of the overall data/MC 
efficiency ratio is simply given by 



R - — Y HijkCijk - YjRijkCi, 



H, 



ijk 



7* 



(29) 



Calculation of the statistical uncertainty on this num- 
ber is slightly non-trivial since we have to take into ac- 
count the fact that the statistical uncertainty on the nor- 
malization bin ratio, ctr,^,, influences the entire row (ij). 
Substituting Rjjk in Eq.|39]with the expression in Eq.lZTl 
and differentiating the resulting expression with respect 
to each of the variables that enter it while using Rij\ - 
1, we obtain 



HijkRjjk + 



(30) 



where the first term reflects the finite size of the signal 
MC sample used in the study and is generally the small- 
est, the second term reflects the statistical uncertainties 
on the number of K^^ 's in bins other than the normaliza- 
tion bin, and the third term, the dominant one, reflects 
the dependence on the statistical precision of Riji of the 
correction factors in each of the bins Rjjk ,ki^\. 

The correction factors are applied to signal MC 
for the decay modes ^ (^/T^' and n^D iD- 
K'-^tt^), which provide K'^ spectra representative of most 
cases of used in BaBar analyses, to determine the 
overall correction factor and its statistical error. The 
above exercise is repeated for several sets of K'^ qual- 
ity cuts, from none to tight, and for three different bin- 
ning approaches. Half of the largest deviation in the K'^ 
correction factors for different quaUty cuts is con- 
sidered to be the systematic uncertainty associated with 
daughter reconstruction efficiency. For these modes 
we are able to determine the ratio of the data and MC 
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1003 daughter reconstruction efficiency to be about99.5% 

1004 and with a statistical error of ~ 0.4% and a systematic 

1005 uncertainty of ~QJ%. 



9. Conclusion 
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In conclusion, we studied the track reconstruction 
efficiency of charged particles in BaBar over a wide 
range of momentum, polar angle, and track separation. 
Our results come from several different control sam- 
ples, which are observed to be self-consistent, and well 
modeled in MC. The overall reconstruction efficiency 
for isolated tracks is found to be consistent with MC 
predictions. We also measured the charge asymmetry 
in the track reconstruction, which was found be con- 
sistent with zero. Any observed difference between 
data and MC in the track reconstruction efficiency could 
be considered as a source of systematic uncertainty in 
all the physics analyses in BaBar. For physics anal- 
yses with low multiplicity and similar topology to the 
Tau31 decays and for charged tracks with momentum 
greater than 180 MeV/c, the results from the Tau31 
study should be used; for other B and D decays, the 
results from the ISR study could be considered. For 
charged tracks with momenta less than 180 MeV/c, an 
additional systematic uncertainty of 1.54% per track 
should be applied. The K'^ reconstruction efficiency in 
data and MC is found to be a function of A'^' momentum, 
polar angle and transverse flight distance, which needs 
to be considered for A"" reconstruction in BaBar. 
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